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C3 , Abstract. We give a functional analytic construction of the fermionic projector on 

a globally hyperbolic Lorentzian manifold of finite lifetime. The integral kernel of the 
£C) ' fermionic projector is represented by a two-point distribution on the manifold. By 

£N| . introducing an ultraviolet regularization, we get to the framework of causal fermion 

systems. The connection to the "negative-energy solutions" of the Dirac equation 
— h , and to the WKB approximation is explained and quantified by a detailed analysis 

• of closed Friedmann-Robertson-Walker universes. 
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1. Introduction 

The fermionic projector was introduced in [7] as an operator which gives a splitting 
of the solution space of the Dirac equation into two subspaces (see also [El Chapter 2] 
and [II])- In a static space-time, these subspaces reduce to the spaces of positive and 
negative energy which are familiar from the usual Dirac sea construction. The signifi- 
cance of the fermionic projector lies in the fact that it can be constructed canonically 
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even in the time-dependent setting. It plays a central role in the fermionic projec- 
tor approach to relativistic quantum field theory (see the review article [10] and the 
references therein). 

So far, the fermionic projector was only constructed perturbatively in a formal 
power expansion in the potentials in the Dirac equation. In the present paper, we give 
a non-perturbative construction of the fermionic projector. To this end, we consider the 
Dirac equation on a globally hyperbolic Lorentzian manifold. For technical simplicity, 
we assume that space-time has finite lifetime. A space-time of infinite lifetime (like 
Minkowski space) can be treated with the same ideas and methods, using the so-called 
mass oscillation property as an additional technical tool. Since the mass oscillation 
property is of independent interest, we decided to work out the case of an infinite 
lifetime in a separate paper [16] . 

In order to explain the basic difficulty which prevented a non-perturbative treatment 
so far, we briefly outline the construction in [7J. Suppose that we consider the Dirac 
equation in Minkowski space (M, (., .)) in a given external potential B, 

(t# + S-m)^ = 0. 

Then the advanced and retarded Green's functions and are solutions of the 
distributional equations 

(i0 + B _ m) s^(x, y) = 5 4 (x - y) = (tf + B - m) s^(x, y) . 

They are uniquely defined by the conditions that the distribution s v (x, .) (and s A (x, .)) 
should be supported in the causal future (respectively past) of x. Taking the difference 
of the advanced and retarded Green's function gives a solution of the homogeneous 
Dirac equation, which we refer to as the causal fundamental solution k m , 

k m (x,y) := — (s v (x,y) - s A (x,y)) . 
We also consider k m as the integral kernel of a corresponding operator 

Jm 

which acts on the wave functions in space-time. Formally, the fermionic projector is 
obtained by taking the absolute value of this operator, 

formally 

Pm ■ — \""m\ s K 1 - 1 ) 



and by forming the combination 

p ( x , V) ■= 5 (Pm(z, y) ~ k m {x, y)) 

(for the rescaling procedure needed to obtain the proper normalization see [II]). The 
basic difficulty is related to the fact that the operator k m acts on the wave functions 
in space-time, which do not form a Hilbert space. More specifically, k m is symmetric 
with respect to the Lorentz invariant inner product on the wave functions 



<ifj\(f)> = / ip(x)<f>(x) (fx (1.2) 
Jm 

(where tp = ip'j is the so-called adjoint spinor). But as (ll.2j) is not positive definite, 
the corresponding function space merely is a Krein space. There is a spectral theorem 
in Krein spaces (see for example [U [20]), but this theorem only applies to so-called 
definitizable operators. The operator k m , however, is not known to be definitizable, 
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making it impossible to apply spectral methods in indefinite inner product spaces. 
The methods in [7\ give a mathematical meaning to the absolute value in (jl.ip in a 
perturbation expansion, leading to the so-called causal perturbation theory. But a 
non-perturbative treatment seemed out of reach. 

We now outline our method for bypassing the above difficulty, again for an external 
potential in Minkowski space. One ingredient is work instead of the space of wave 
functions with the solution space of the Dirac equation. This solution space has a 
natural Lorentz invariant scalar product 

(</#) := / (? 7 V)(*,aO d 3 x , (1.3) 

giving rise to a Hilbert space 9i m . Our main observation is that the operator k m 
relates the scalar product (|1.3|) to the space-time inner product (|1.2|) by 

(^\k m <p) = <^\4» (1.4) 

(valid if ift is a solution of the Dirac equation; see Proposition 13. II below) . On the other 
hand, we can express the bilinear form <.|.> in terms of the scalar product using a 
signature operator S, 

<i/}\<j» = (if)\S<l>) (1.5) 

(valid if ip and (ft are solutions of the Dirac equation; see equation (|3.3p below). 
The operator S will turn out to be a bounded symmetric operator on the Hilbert 
space (JK m , (■{■))■ Comparing (jl.4p with (|1.5p . we find that on solutions of the Dirac 
equation, the operator k m can be identified with the operator S. This makes it possible 
to use spectral theory in Hilbert spaces to define the absolute value in (jl.ip . 

In Section [3j we will make this construction mathematically precise in the setting of a 
globally hyperbolic space-time of finite lifetime. We point out that all our constructions 
are manifestly covariant. They do not depend on the choice of a foliation of the 
manifold. It makes no difference whether the Cauchy surfaces are compact or non- 
compact. We do not need to make any assumptions on the asymptotic behavior of the 
metric at infinity. 

In Section HI it is explained how the fermionic projector gives rise to examples of 
causal fermion systems as defined in [131 Section 1]. 

Our construction of the fermionic projector gives a splitting of the solution space 
of the Dirac equation into two subspaces. For the physical interpretation, it is im- 
portant to understand how these subspaces relate to the usual concept of solutions 
of positive and negative energy. To this end, we analyze the fermionic projector in a 
closed Friedmann-Robertson- Walker universe. This has the advantage that the Dirac 
equation reduces to an ODE in time, which can be analyzed in detail. In particular, 
the concept of "solutions of negative energy" (which for clarity we mostly refer to as 
"solutions of negative frequency") can be made precise by a specific WKB approx- 
imation as worked out in [15] . In Section \5\ it is shown that our definition of the 
fermionic projector agrees with the concept of "all solutions of negative frequency," 
provided that the metric is "nearly constant" on the Compton scale as quantified in 
Theorem 15.11 and Theorem 15.21 It is remarkable that, in contrast to a Gronwall es- 
timate, our error estimates do not involve a time integral of the error term. This 
means that small local errors of the WKB approximation do not "add up" to give a 
big error after a long time. Moreover, our estimates also apply near the big bang and 
big crunch singularities. Keeping these facts in mind, our estimates show that for our 
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physical universe, the fermionic projector coincides with very high precision with the 
usual concept of the Dirac sea being composed of all negative-frequency solutions of 
the Dirac equation. This gives a rigorous justification of the physical concepts behind 
the fermionic projector approach. 

In Section [6l we analyze what happens if the metric changes substantially on the 
Compton scale. To this end, we consider a closed Friedmann-Robertson- Walker uni- 
verse with a scale function R(r) being piecewise constant. Similar to the situation for 
Klein's paradox, at the times when R is discontinuous, the frequencies of the solutions 
change. As a consequence, the concept of positive or negative frequency becomes mean- 
ingless. In this situation, our constructions still apply, giving a well-defined fermionic 
projector. This fermionic projector consists of a mixture of positive and negative fre- 
quencies. Moreover, as we explain in an explicit example where S = 0, the fermionic 
projector may depend sensitively on the detailed geometry of space-time. 

2. Preliminaries 

Let (M,g) be a smooth, four-dimensional, globally hyperbolic Lorentzian manifold. 

For the signature of the metric we use the convention (H ). As proven in [3], M 

admits a smooth foliation (N t )t£R by Cauchy hypersurfaces. Thus M is topologically 
the product of R with a three-dimensional manifold. This implies that M is spin 
(for details see [21 EI])- We let SM be the spinor bundle on M and denote the 
smooth sections of the spinor bundle by C°°(M, SM). Similarly, C^°(M, SM) denotes 
the smooth sections with compact support. The fibres S X M are endowed with an 
inner product of signature (2,2), which we denote by -<.\.y x . The Lorentzian metric 
induces a Levi-Civita connection and a spin connection, which we both denote by V. 
Every vector of the tangent space acts on the corresponding spinor space by Clifford 
multiplication. Clifford multiplication is related to the Lorentzian metric via the anti- 
commutation relations. Denoting the mapping from the tangent space to the linear 
operators on the spinor space by 7, we thus have 

7 : T X M->L(S X M) with 7(11) j{v) + i(v) 7(1*) = 2g(u, v) ls x (M) ■ 

We also write Clifford multiplication in components with the Dirac matrices 7^ and 
use the short notation with the Feynman dagger, j(u) = vP^j = t/l. The connections, 
inner products and Clifford multiplication satisfy Leibniz rules and compatibility con- 
ditions; we refer to [21 EI] for details. Combining the spin connection with Clifford 
multiplication gives the geometric Dirac operator D = ij^Vj. In order to include 
the situation when an external potential is present, we add a multiplication opera- 
tor "B(x) G h(S x M), which we assume to be smooth and symmetric with respect to 
the spin scalar product, 

B G C°°(M,L(SM)) with -<S0|^^x = ^\'Bipy x G S X M . (2.1) 

We then introduce the Dirac operator by 

V := ij j Vj + S : C°°(M, SM) -> C°°(M, SM) . (2.2) 

For a given real parameter m G R (the "rest mass"), the Dirac equation reads 

{V-m)i; m = 0. (2.3) 

For clarity, solutions of the Dirac equation always carry a subscript m. 
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In the Cauchy problem, one seeks for a solution of the Dirac equation with initial 
data ipN prescribed on a given Cauchy surface N. Thus in the smooth setting, 

(D-m)^ m = 0, ^\ N = ^ N eC°°(N,SM). (2.4) 

This Cauchy problem has a unique solution tp m £ C°°(M,SM). This can be seen 
either by considering energy estimates for symmetric hyperbolic systems (see for ex- 
ample [12]) or alternatively by constructing the Green's kernel (see for example [T]). 
These methods also show that the Dirac equation is causal, meaning that the solution 
of the Cauchy problem only depends on the initial data in the causal past or future. 
In particular, if ipjy has compact support, the solution tp m will also have compact sup- 
port on any other Cauchy hypersurface. This leads us to consider solutions ip m in the 
class C^?(M, SM) of smooth sections with spatially compact support. On solutions in 
this class, one introduces the scalar product (.|.)jv bjQ 

(^m\4>m)N = 27T / ^V>m|?%n^z d/Mr(x) , (2.5) 



JN 

where f denotes Clifford multiplication by the future-directed normal v (we always 
adopt the convention that the inner product -<.\^.y x is positive definite). This scalar 
product does not depend on the choice of the Cauchy surface N. To see this, we let N' 
be another Cauchy surface and Q the space-time region enclosed by N and N'. Using 
the symmetry property in (12. ip together with (12 . 2 [) and (12.3|) . we obtain 

iVi^VmiyVm^ic = -<(—iVj)lpm\7 i <l>m>~x + -<^m\{i'f^j)^m>~x , n 

(2.DJ 

= <Vtp m \(j) m >- x + ^1prn\V(f) m >- x = , 

showing that the vector field -Kip m \'y- 1 ] <j>m>~x is divergence- free ("current conservation"). 
Integrating over Q and applying the Gaufi divergence theorem, we find that (i/j m \<fim)N = 
("0m I ^m) N' ■ I* 1 view of the independence of the choice of the Cauchy surface, we simply 
denote the scalar product (|2.5p by (.|.). Forming the completion, we obtain the Hilbert 
space (IK m ,(.|.)). It consists of all weak solutions of the Dirac equation (|2.3p which 
are square integrable over any Cauchy surface. 

The retarded and advanced Green's operators s!^ and are mappings (see for 
example [I]) 

s^s^ : C%°(M,SM) -> C™(M,SM) . 

They satisfy the defining equation of the Green's operator 

{V-m){ S ^<p) =4>. (2.7) 

Moreover, they are uniquely determined by the condition that the support of s^cj) 
(or s^cj)) lies in the future (respectively the past) of supp <j). The causal fundamental 
solution k m is introduced by 

km := ^~ {s m ~ stk) ■ C^{M, SM) -> C™(M, SM) n "K m . (2.8) 

It gives rise to an explicit solution of the Cauchy problem, as we recall in the next 
lemma. We only sketch the proof, because in Lemma [3.91 an independent proof will be 
given. 



^The factor 2-k might seem unconventional. This convention was first adopted in [13]. It will 
simplify many formulas in this paper. 
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Lemma 2.1. The solution of the Cauchy problem (|2.4p has the representation 

ip m (x) = 2ir k m (x,y)i/ip N (y) dfi N (y) , 
Jn 

where k m (x,y) is the integral kernel of the operator k m , i.e. 

(k m <j))(x) = / k m (x,y)cf)(y) dfj, M (y) ■ (2.9) 
Jm 

Sketch of the Proof. For the proof that k m can be represented with an integral ker- 
nel (|2,9p and for analytic details on k m (x,y) we refer to pp. In order to prove (|2.4|) . it 
suffices to consider a point x in the future of N, in which case fj2.4j) simplifies in view 
of (USD to 

i/) m (x)=i s^ n (x,y)i/'(y)ip N (y) dfi N (y) . 
Jn 

This identity is derived as follows: We let r/ 6 C°°(M) be a function which is identically 
equal to one at x and on N, but such that the function r\i\) m has compact support (for 
example, in a foliation (N t )t<^R one can take rj = x(t) with \ £ Co° (!&))• Then, 
using (HZD, 

*Ms) = Mm)(z) ( = } s*((0) - m)(# TO )) = shW&w)*™)) > ( 2 - 10 ) 

where we used (|2.7p and the fact that V>m is a solution of the Dirac equation. In (*) 
we used the identity 

Tp = s A ((D-m)Tp) for ip G C^°(M, S'M) , 

which follows from the uniqueness of the solution of the Cauchy problem, noting that 
the function tp — s A ((T> — m)ip) satisfies the Dirac equation and vanishes in the past 
of the support of ip. To conclude the proof, as the function r\ in fj2. 10f> we choose a 
sequence rji which converges in the distributional sense to the function which in the 
future and past of N is equal to one and zero, respectively. □ 



3. Functional Analytic Construction of the Fermionic Projector 

3.1. The Space-Time Inner Product as a Dual Pairing. On the Dirac wave 
functions, one can introduce the Lorentz invariant inner product 

<V#> := / ^4>dfi M ■ (3.1) 
Jm 

In order to ensure that the space-time integral is finite, we assume that one factor has 
compact support. In particular, we can regard <.|.> as the dual pairing 

<.|.> : IK m x C$°(M,SM) -> C . 

The next proposition shows that the causal fundamental solution is the signature 
operator of this dual pairing. 

Proposition 3.1. For any i(j m £ "K m and <p £ C°°{M,SM), 

(ipm \km4>) = <V>m|</» • 
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Proof. We choose Cauchy surfaces N + and iV_ lying in the future and past of supp cp, 
respectively. Let fi be the space-time region between these two Cauchy surfaces, i.e. 
dQ = N + U 2V_. Then, according to $ZM . 

(tpm \km4>) = (tp m \ k m (j)) N+ = — (lp m \ (p) N+ 

i r 

= i \ Vj^tpm | j j s^(j)y x dfi{x) , 

where in the last line we applied the GauB divergence theorem and used (|2,5p . Using 
that tp m satisfies the Dirac equation, a calculation similar to (j2.6|) yields 

(lpm\krn4 t )= / \ (P ~ m) S^y x dfi(x) = / -<tp m \^»~x dfJ,(x) . 

Jn Jn 
As <j> is supported in Q, we can extend the last integration to all of M, giving the 
result. □ 



Corollary 3.2. The operator k rn , (|2.8p . is symmetric with respect to the inner prod- 
uct (IP). 



Proof. Using Proposition 13.11 we obtain for all 0, ifi S Cq*(M,SM), 

<k m (j) | -0> = (^m0 | k m i>) = <(f) | fc m ^> , 

concluding the proof. □ 

3.2. Space-Times of Finite Lifetime. For the construction of the fermionic pro- 
jector, we need to assume that space-time has the following property. 

Definition 3.3. A globally hyperbolic manifold (M,g) is said to be m-finite if there is 
a constant c > such that for all 4> m ,ip m £ 5{ m nC^(M, SM), the function «j) m \ipm>~x 
is integrable on M and 

\<4>m\^m>\ < c \\<f> m \\ \\ip m \\ (3.2) 
(where \\.\\ = (.|.)2 is the norm on Jt m ). 

Before going on, let us briefly discuss which manifolds are m-finite. 

Definition 3.4. A globally hyperbolic manifold (M, g) has finite lifetime if it admits 
a foliation (^Vt)te(*0)*x) Cauchy surfaces with a bounded time function t such that the 
function (v, t? t ) is bounded on M (where v denotes the future- directed normal on N t 
and {u,d t ) = g(u,d t )). 

Proposition 3.5. Every globally hyperbolic manifold of finite lifetime is m-finite. 

Proof. Let ip m £ C™(M, SM) be a solution of the Dirac equation (|2.3p . Applying 
Fubini's theorem and decomposing the volume measure, we obtain 

<1pm\^m> = [ -<1p m \lp m y{x) dfi M {x) = [ [ <^m\^m>- (v, d t )dtdfl Nt 
JM Jt JN t 

and thus ^ 

\<1pm\lpm>\ < SUp(u, d t ) / dt \^1p m \^ m ^\d^ Nt . 
M Jt JN t 
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Estimating the spatial integral by 

HVVilVVnH dn Nt < / -<1p m \j">Pm>- dfiNt = (VVnlVw) > 



we conclude that 



\<1pm\4>m>\ < {h -t ) 8Up{v,d t ) 

M 



Polarization and a density argument give the result. □ 

Proposition 3.6. On a globally hyperbolic manifold of finite lifetime, there is a con- 
stant C < oo such that the arc length of every timelike geodesic is at most C . 

Proof. Let 7 be a timelike geodesic. Possibly after extending it, we can parametrize 
it by the time function t £ (to,ti) of our foliation. Then the vector field 7 — dt is 
tangential to Nt- Hence we can estimate the length of the geodesic by 

L{l)= r \f{i,l)dt< f 1 y / (^,u)(u^)dt = [\u,d t )dt < (ti-to)sup(i/,^) . 



'to Jt Jt M 

This concludes the proof. □ 

We do not know whether an upper bound on the length of timelike geodesies already 
implies that the space-time has finite lifetime in the sense of Definition 13.41 Moreover, 
it is not known if every m-finite manifold necessarily has finite lifetime. Unfortunately, 
entering the study of these open questions goes beyond the scope of the present paper. 

3.3. Definition of the Fermionic Projector. Let us assume that (M,g) is m-finite. 
Then the space-time inner product can be extended by continuity to a bilinear form 

< . I . > : 3~C m x !J~C m — > C . 

Moreover, applying the Riesz representation theorem, we can uniquely represent this 
inner product with a signature operator S, 

S : 'K ra ^r c K m with «j) m \^m,> = (4>m I S ipm) . (3.3) 



The operator S is obviously symmetric. Moreover, it is bounded according to (|3.2p . We 
conclude that it is self-adjoint. The spectral theorem gives the spectral decomposition 



S= XdE x , 

where E\ is the spectral measure (see for example |24j). The spectral measure gives 
rise to the spectral calculus 



f(S) = [ /(A) dE x , 



la(S) 

where / is a Borel function. 

The spectral calculus for the operator S is very useful because it gives rise to a 
corresponding spectral calculus for the operator k m , as we now explain. Multiplying 
k m from the left by f(S) with a bounded function / G C°(a(S),M.) gives an operator 

f(S)k m : C yD (M, SM) — > "K m . 
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This operator is again symmetric with respect to <.|.>, because for any <fi,ip £ 
C^(M,SM), 

<f(S)k m ^)\ip> = (f(S)k m (f>\k m ip) 

= (km<t>\f(S)k m tp) = <(j>\f(S)k m ip> , 

where in the first and last equality we applied Proposition 13. 11 In order to make sense 
of products of such operators, we can consider the inner product <f(S)k m (p\g(S)k m 'ip> 
with /, g € C°((t(5),R). Combining (I3.3P with the spectral calculus for S and Propo- 
sition [3Jl we obtain 



<f(S)k m (j)\g(S)k m 'ip> = (f(S)k m (j)\S g(S)k m ip) 

= (k m 4> \ (fg) (S) Sk m ij)) = \ (fg) {S)Sk m ^>. (3 ' 5) 
In view of (|3.4|) . this identity can be written in the suggestive form 

(f(S)k m ) (g(S)k m ) for ^ ally (fg)(S) S k m . (3.6) 

Note that this last equation makes no direct mathematical sense because the image 
of the operator g(S)k m does not lie in the domain of k m , making it impossible to 
take the product. However, with fj3.4|) and ()3.5[) we have given this product a precise 
mathematical meaning. 

We now use this procedure to construct the fermionic projector. 

Definition 3.7. Assume that the globally hyperbolic manifold (M,g) is m-finite (see 
Definition\EM. Then the operators P± : C$°(M,SM) -> "K m are defined by 

P+ = X(o,oo)(S) k m and P- = -X(-oo,o)(S) k m (3.7) 

(where \ denotes the characteristic function). The fermionic operator P is defined 
by P = P_. 

Proposition 3.8. For all 4>, ifi £ Cq°(M, SM), the operators P± have the following 
properties: 

<P± <j> | ip> = <(j> | P± tp> (symmetry) (3.8) 

<P + 4> | P- 4>> = (orthogonality) (3.9) 

<P± 4> I Pt i» = <4> I \S\P± 4» (normalization). (3.10) 

Moreover, the image of P± is the positive respectively negative spectral subspace of S, 
meaning that 



P+(C °°(M, SM)) = P (0 ,oo)(^ m ) , P-(C °°(M, SM)) = E^flpin) 



Proof. This follows immediately from (|3.4p . (|3.5p and the functional calculus for self- 
adjoint operators in Hilbert spaces. □ 

We finally explain the normalization property (|3.10p . We first point out that, due to 
the factor S on the right of (|3.10p . the fermionic operator is not a projection operator. 
The projection property could have been arranged by modifying (|3.T[) to 

P = ~X(-oo,o)(S) 1ST 1 k m . 

However, we prefer the definition (|3.7p and the normalization (I3.10P for the following 
reason. In the perturbative construction in Minkowski space (HE], we worked with a 
^-normalization in the mass parameter (for details see eqns (3.19)-(3.21)] or |llj). 
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This ^-normalization was the original motivation for the nomenclature "fermionic pro- 
jector." Clearly, such a 5- normalization in the mass parameter cannot be used in a 
space-time of finite lifetime where all appearing space-time integrals are finite. But 
the (5-normalization in the mass parameter will again be used in the non-perturbative 
construction on a globally hyperbolic space-time of infinite lifetime |16] , In order to get 
agreement with the formulas in |16| . we must work with the normalization as in (|3.7p . 
In order to avoid confusion, we here call the operator P the fermionic operator. This 
also harmonizes with the notions in the framework of causal fermion systems, as will 
be explained in Section H] below. 

3.4. Explicit Formulas in a Foliation. It is instructive to supplement the previous 
abstract constructions by explicit formulas in a foliation. We always work with the 
following particularly convenient class of foliations. As shown in [4, 22], there are folia- 
tions (iVj)tgR by Cauchy surfaces where the gradient of the time function is orthogonal 
to the leaves and the lapse function is bounded, i.e. 

g = (3 2 dt 2 - g Nt with < (3 < 1 , (3.11) 

where gN t is the induced Riemannian metric on Nt, and the lapse function (3 is a smooth 
function on M. We remark that in space-times of finite life time (see Definition [33]), the 
time parameter t could be chosen on a bounded interval. In this case, for convenience 
we prefer to parametrize t on all of R, such that lim /3 = 0. We denote space-time 

t— >±oo 

points by (t,x) with t £ M and x E N t . Moreover, we denote the scalar product (|2.5|) 
for N = N t by (.|.)t, and the corresponding Hubert space by (!K t ,N t ). Solving the 
Cauchy problem with initial data on Nt and evaluating the solution at another time t' 
gives rise to a unitary time evolution operator 

: JCt -> "K t > ■ 

Clearly, the unitary time evolution operators are a representation of the group (R, +). 
The time evolution also gives rise to the unitary mapping 

which allows us to canonically identify each Hilbert space (5Cj, (.|.)t) with (IK m , (.|.)). 
We denote the restriction of a smooth wave function ip £ C oc (M,SM) to the hyper- 
surface N t by tpu. 

Lemma 3.9. For every (f> £ Cfi°(M,SM), 

(s^)(t,x) = -i ! (U^^f<f, w )){x) dt' (3.12) 

(k m <j>){t,x) = ^j°° (U^\^ tl )){x) dt' . (3.13) 

Proof. The Dirac operator can be written as 

V = /TV (id t - H t ) , 

where Ht is a purely spatial operator acting on "Kt (the "Hamiltonian" ) . We apply 
the Dirac operator to the right side of (|3.12j) . which we denote by F(t,x). As the 
integrand in (|3.12p is a solution of the Dirac equation, only the derivative of the limit 
of integration needs to be taken into account, 

(V - m)F(t, x) = (/TV(i, x)) (U^iP^t)) (x) . 
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Using that £/*'* is the identity, we conclude that 

(V - m)F(t,x) = <f>(t,x) . 

Hence F(t,x) satisfies the defining equation of the Green's operator (|2.7p . Moreover, 
it is obvious that F(t, x) vanishes if t is in the past of the support of 0. The unique 
solution of the Cauchy problem gives the result. 

Repeating the above argument for the advanced Green's operator gives 

(sU)(t, x ) =i j™ (^'*'(/¥<V)) (*) dt> ■ 

We finally apply ([278]) to obtain ||515)> . □ 

For what follows, it is useful to identify !K m with the Hilbert space "Kt for some 
fixed time to. The formulas of the previous lemma are then rewritten by multiplying 
with the time evolution operator. For example, 

1 f°° 

k m 4> = C/ t0 '*(/3^) |t dt : C °°(M, SM) -> X t0 . (3.14) 

Lemma 3.10. The operator S as given by fj3.3j) has the representation 

to ■ 



1 f°° 

S=— U tQ t (Pf), t U^ dt : X to ->• %, 

27r i-oo 



Proof. Rewriting the space-time integral in (|3.ip with Fubini's theorem and using the 
identity ^ = 1, we obtain 

«t>m\^ m > = J {^j -«Pm\ip m y (t,x) P(t,x) dfj, Nt (x)J dt 
1 f°° 

7j~ Wm | Wf)\t^m)t dt 



2tt 
1 

2^ 



-oo 
oo 



-oo 

Comparing with (|3.3p gives the result. □ 



Iterating (|3.13|) . we can make the following formal calculation, 

(k m k m <p) {t0 = — dt dt'U^^f)^'^)^ 



— oo J — oo 

OO /'OO 



1 />oo /'OO 

^ / ^ / dt' u^iMh ut,t0 ^(WV 



OO J — OO 



where in the second line we used the group property of the time evolution operator. 
Comparing with (|3.14p . we obtain the simple relation 

t , formally 

This is precisely the relation (|3.6p in the special case f,g= 1. Iteration gives similar 
formal expressions for polynomials of k m , from which (j3.6f) can be obtained formally 
by approximation. Although the last arguments are only formal, they explain how 
the functional calculus (|3.6p comes about. In order to give this functional calculus 
a mathematical meaning, one needs to evaluate weakly as is made precise by ([37 
and (13751). 
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3.5. Representation as a Distribution. We now represent the fermionic operator 
by a two-point distribution on M . 

Theorem 3.11. There is a unique distribution V G D'(MxM) such that for all(ft,ip G 
C$°(M,SM), 

«P\PiP> = V{(j)®^) . 
Proof. According to Definition 13.71 

\\P(<t>)\\ = IIX(-oo,o)(5 , )M0)ll- 
Since the norm of the operator X(-oo,o) {&) ls bounded by one, we conclude that 

\\P(m < \\km<t>\\ ■ (3.15) 
Using ()3.13|) . we can estimate the last norm in a foliation (|3.1ip by 
1 f°° ,, 

\\k m <t>\\<— \W<p\t> L dt'<Cm&x\<f>\ with C = C(supp0). (3.16) 
2vr J_ 00 M 

Suppose that <p n G Cq°(M, SM) is a sequence which converges to zero in the sense 
that there is a compact set K C M with supp (j) n C K for all n and that 

Emax \d a d>ri\ n ~ > °°> for some N > . 
|o|<jv 1 rn| - 

Then the estimates f)3. 15|) and (|3,16p show that P{4> n ) — > in 5C m . Taking the scalar 
product with the vector k m ^ip with ip G C^°, we conclude that 

= lim (P((pn) I k m i/>if) = lim <P(<^„) | i/>ip> = lim / -<P(0 n )|^>- x d^ M (^) • 

Since the inner product is positive definite and ip is arbitrary, we conclude 

that P(4> n ) converges to zero in L\ oc (M). Hence P(4> n ) — > as a distribution. 

The result now follows from the Schwartz kernel theorem (see [181 Theorem 5.2.1], 
keeping in mind that this theorem applies just as well to bundle- valued distributions 
on a manifold simply by working with the components in local coordinates and a local 
trivialization). □ 

In order to get the connection to [8j, it is convenient to use the standard notation 
with an integral kernel P(x,y), 

«t>\Pip> = -<H x ) I p( x , v) ^(y)>~x dn M {x) dfi M (y) 

J J MxM 



(Pi>)(x) = / P(x,y)i;(y)d m (y) . 

JM 

In view of Proposition 13. 8[ we know that last integral is not only a distribution, but a 
function which is square integrable over every Cauchy surface. Moreover, the symmetry 
of P, (I3.8p . implies that 

P(x,y)*=P(y,x), 
where the star denotes the adjoint with respect to the spin scalar product. 
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4. Connection to the Framework of Causal Fermion Systems 

We now explain the relation to the framework of causal fermion systems as intro- 
duced in [13] (see also [E]). In order to get into this framework, we need to introduce 
an ultraviolet regularization. This is done most conveniently with so-called regular- 
ization operators. 

Definition 4.1. A family (9^ e ) e >o of bounded linear operators on ?{ m are called reg- 
ularization operators if they have the following properties: 

(i) Solutions of the Dirac equation are mapped to continuous solutions, 

V\ £ : Jim -> C°(M, SM) n W m 

(ii) For every e > and x G M, there is a constant c > such that 

\\(^ m )(x)\\ < c\\^ m \\ V*p m ex m . (4.1) 

(where the norm on the left is any norm on S X M). 

(iii) In the limit £ \ 0, the regularization operators go over to the identity with 
strong convergence of$K e and9\*, i.e. 

JH.Vm, K^m ^ ?Pm »« X m V 1p m G !K m . (4.2) 

There are many possibilities to choose regularization operators. As a typical exam- 
ple, one can choose finite-dimensional subspaces < K^ > C C££(M, SM) n !K m which are 
an exhaustion of "K m in the sense that Ji^ c jtf 1 ) C ••• and 3i m = iyK( n ). Set- 
ting n(e) = max([0, 1/e] n N), we can introduce the operators 5H £ as the orthogonal 
projection operators to "K^ n ^ . An alternative method is to choose a Cauchy hyper- 
surface N, to mollify the restriction VViliV to the Cauchy surface on the length scale e, 
and to define 9\ e ip m as the solution of the Cauchy problem for the mollified initial 
data. 

Given regularization operators 9l e , for any e > we introduce the particle space 
(^particle, (-lOrKpartidJ as the Hilbert space 

^particle = ker (*H e X(-QO,0) (S)) ) (■ I -^particle = (' I ') l^particlo x:K particlc " 

Next, for any x G M we consider the bilinear form 

U . ^H^particlc ^ -^particle 

This bilinear form is bounded in view of (I4.ip . The local correlation operator F £ (x) is 
defined as the signature operator of this bilinear form, i.e. 

b(tp m , (j> m ) = {ipm \ F £ (x) <j) m ) J{ particlc for all 1p m , (fi m 

Taking into account that the spin scalar operator has signature (2,2), the local cor- 
relation operator is a symmetric operator in L(5C particlc ) of rank at most four, which 
has at most two positive and at most two negative eigenvalues. Finally, we introduce 
the universal measure p = F £ dpM as the push-forward of the volume measure on M 
under the mapping F £ (thus p(fl) := /iM((^ e ) _1 (^))- Omitting the subscript "parti- 
cle", we thus obtain a causal fermion system of spin dimension two as defined in [13l 
Section 1.2]: 

Definition 4.2. Given a complex Hilbert space ("K, (.|.)jc) "particle space",) 
and a parameter n G N (the "spin dimension" ), we let 3~ C L(!K) be the set of all 
self-adjoint operators on "K of finite rank, which (counting with multiplicities) have 
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at most n positive and at most n negative eigenvalues. On 3~ we are given a positive 
measure p (defined on a a -algebra of subsets of J), the so-called universal measure. 
We refer to (!K, 3~, p) as a causal fermion system in the particle representation. 

The formulation as a causal fermion system gives contact to a general mathemati- 
cal framework in which there are many inherent analytic and geometric structures 
(see [9j Q2] )■ In particular, the differential geometric objects of spin geometry have 
a canonical generalization to the regularized theory. Namely starting from a causal 
fermion system (IK, 3~, p) one defines space-time as the support of the universal mea- 
sure, M := suppp. Note that with this definition, the space-time points x,y £ M are 
operators on K (thinking of our above construction of the causal fermion system, this 
means that we identify a space-time point x with its local correlation operator F £ (x)). 
On M, we consider the topology induced by 9" C L(K). The causal structure is 
encoded in the spectrum of the operator products xy: 

Definition 4.3. For any x,y £ 3~, the product xy is an operator of rank at most 2n. 
We denote its non-trivial eigenvalues by , . . . , A^ (where we count with algebraic 
multiplicities). The points x and y are called timelike separated if the Aj are all real. 
They are said to be spacelike separated if all the X* y are complex and have the same 
absolute value. In all other cases, the points x and y are said to be lightlike separated. 

Next, we define the spin space S x by S x = x(K) C K endowed with the inner prod- 
uct -<-\->- x := —{.\x.)j{. The kernel of the fermionic operator with regularization is 
introduced by 

P e (x,y) = ir x y : S y -> S x , (4.3) 
where tt x is the orthogonal projection to S x in K. Connection and curvature can be 
defined as in |12t Section 3]. We remark for clarity that the Dirac equation and the 
bosonic field equations (like the Maxwell or Einstein equations) cannot be formulated 
intrinsically in a causal fermion system. Instead, as the main analytic structure one 
has the causal action principle. 

We conclude this section by deriving more explicit formulas for the local correlation 
operators. Moreover, we compute the regularized fermionic operator and compare it 
to the unregularized fermionic operator of Definition 13.71 To this end, for any x G M 
we define the evaluation map e e x by 

e e x : K m -> S X M , e £ x i> m = (<H e X(-oo,o)(S) 1> m ){x) . (4.4) 
We denote its adjoint by l x , 

4 := (4)* : S X M -> <K m . 
Multiplying i e x by e e x gives us back the local correlation operator F e (x) (extended by 
zero to the orthogonal complement of IKparticio), 

F £ (x) = -4 e £ x : K m -> K m . (4.5) 

Let us compute the adjoint of the evaluation map. For any \ £ S X M and ip m G !H m , 
we have according to (|4.4p 

((4)*X|^m) = <X\^eX{-oofi){ S )^m^x = <SxX\%>X(-oc,0)(S) V>m> , 

where 5 X is the 5-distribution supported at x (thus in local coordinates, 5 x (y) = 
| det g(x)\~z 5 A (x — y)). Applying Proposition 13.11 gives 

((4)*XlV'm) = (k m S x x\^eX(-oo,0)(S)lpm) = (X(-oo,0) ( S ) k m 0~x X I i>m) , 
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and thus 

4 = (4)* = X(-oo,o)(S) K k m S x . (4.6) 

Combining this relation with (|4.4p and (|4.5|) . the local correlation operator takes the 
more explicit form 

F £ (x) = -i%e% = -X(-oo,o){S)^Kk m 6 x m £ X(-oo,o)(S) . 

We next introduce the kernel of the regularized fermionic operator by 

P £ (x,y) = -e £ x L £ y . (4.7) 

After suitably identifying the spinor spaces S X M and S y M with the corresponding 
spin spaces S x and S y , this definition indeed agrees with the abstract definition (|4.3p 
(for details see [121 Section 4.1]). Even without going through the details of this 
identification, the definition (|4.7p can be understood immediately by computing the 
eigenvalues of the closed chain. Starting from the definition (|4.3p . the corresponding 
closed chain is given by A e xy := P £ (x,y) P £ (y,x) = ir x yx7Ty. Keeping in mind that 
in (|4.3p the space-time points are identified with the corresponding local correlation 
matrices, this means that the spectrum of the closed chain is the same as that of the 
product F{y)F{x) (except possibly for irrelevant zeros in the spectrum). Taking the 
alternative definition (|4.7p as the starting point, the closed chain is given by 

A £ {x,y) = (4 L £ y ) (e £ y Q . 

Since a cyclic commutation of the operators has no influence on the eigenvalues, we 
conclude that the closed chain is isospectral to the operator 

44 L x e x = F (y) F ( x ) » 

giving agreement with the abstract definition ()4.3p . 

The corresponding regularized fermionic operator is defined by 

(P e (</>))(*)= / P £ (x,y)<p(y)d m (y). 
JM 

Using (|4.7p together with (|4.6p and (|4.4p . this operator can be written as 

P £ = -m £ X(-oo,o)( S ) K k m ■ C °°(M, SM) C7°(M, SM) n K m . (4.8) 

The next proposition shows that if the regularization is removed, the operator P £ 
converges weakly to P. 

Proposition 4.4. For every <j>,ip G Cfi°(M,SM), 

<0|P £ ^> <(/)\Pij> . 
Proof. Applying Proposition 13. II and (|4.8p . we get 

<0ip £ v» = -(^i^x(-oc,o)^)^:^) = -(^ wix(-oo,o)(5)j«: k*i>)- 

Now use that the operators R* converge strongly according to (14. 2p . □ 
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5. Example: A Closed Friedmann-Robertson- Walker Universe 



We now want to complement the abstract construction of the fermionic projector by 
a detailed analysis in a closed Friedmann-Robertson- Walker space-time. In so-called 
conformal coordinates, the line element reads 



ds 2 



R(t) 2 (dr 2 - d X 2 - sin(x) 2 (d$ 2 + sin 2 dip 2 )) . 



(5.1) 



Here r G (0, tt) is a time coordinate, ip G [0, 2tt) and $ G (0, tt) are angular coordinates, 
and x ^ (0) 7r ) is a radial coordinate. The scale function R(t) should have the following 
properties. We assume that r = and r = tt are the big bang and big crunch 
singularities, respectively. This implies that 

R(0)=0 = R{tt) and i?| (0j7r) >0. 

Moreover, we assume that R is a piecewise monotone C 2 -function (i.e., the inter- 
val (0, tt) can be divided into a finite number of subintervals on which R is monotone). 
It is convenient to write the scale function as 



R(t) = i? max g(r) with 



Ry, 



maxi? 

(0,77) 



(5.2) 



A special case is the dust matter model R(t) = -R max (1 — cos(t)) (see [TTJ, Section 5.3]). 

The spatial dependence of the Dirac equation can be separated by eigenfunctions of 
the Dirac operator on S 3 corresponding to the eigenvalues A G {±§, ±§, • • •} (for details 
see [E]). After this separation, the time evolution operator U T,T ° G C 1 ((0, 7r), U(C 2 )) 
of the Dirac equation is given as the solution of the initial value problem 



mR(r) 



U T0 ' T0 = 1 2 . 
According to Definition 13.71 and ([3 

K(— oo 

i r 

2^ Jo 





-1 



A 



p 



as well as (|3.14p . we have 

0)(S) k m 

(U T > T o\* 





-1 



(r) R(t) dr , 



(5.3) 
(5.4) 

(5.5) 
(5.6) 



where G C§°((0, tt), C 2 ). 

In the subsequent estimates, we shall work with the WKB approximation introduced 
as follows (for more details see [IS])- We first define V(r) as a unitary matrix which 
diagonalizes the coefficient matrix in (|5,3p . i.e. 



V 



Rm 
-A 



-A 
-Rm 



-i 



/ 





-1 



where 



Next, we set 



/(r) := y/X 2 + m 2 R(T) 2 



T,T0 

WKB 



V(t) 



-i 



exp 



V 



exp 



/ 



TO 



V(t ) 



(5.7) 
(5.8) 

(5.9) 



Note that for all t, tq G (0, 7t), the matrices U T ' T °, V(r) and (7^ B are unitary. 
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Applying Lemma 13.101 the signature operator S as denned by (|3.3p takes the form 
S = £v%> r (l ^U^R{r)dr. (5.10) 

Replacing the time evolution by the WKB approximation, we obtain the signature 
operator 

SWKB = jTc^ffl(j ^U^ B R(r)dr. (5.11) 



In analogy to (|5.5p and (|5.6p . we introduce the fermionic projector in the WKB ap- 
proximation by 

AvKB = — X(-oo,0)(5'wKb) fewKB (5-12) 

kwMt) = ± £ (U%£ B y (I ^<t>(r)R(r)dr. (5.13) 

In the following two theorems, we specify under which conditions and in which sense 
the fermionic projector is well-approximated by WKB wave functions. We first state 
the theorems and discuss them afterwards. 

Theorem 5.1. For given tq £ (0,vr) and a given function g, the function Pwkb a s 
defined by (|5.12p can be represented for any values of the parameters A, m and R ma x 
by 

1_ /'%-, .-, /" , ",r, x \ ( 1 

2vr 




1 + 

Theorem 5.2. For any constant k > 0, i/iere is a constant c (only depending on k, 
tq and the function g ), such that for all m and R max with mi? max > 1 the following 
statement holds: For every A in the range 

\X\<kmR max (5.14) 

and every <f> £ C^°((0, tt), C 2 ), we have the estimate 

\\(P - Pwkb) (tfOH <c(mi? max )-l R max I \\<j>{j)\\ dr . (5.15) 



Comparing the exponential factors in (|5.9p with those in Theorem 15. 1| one sees that 
Pwkb only involves the factor exp(i J /), whereas the factor exp(— i J f) in ()5.9[) has 
disappeared. In this sense, our formula of Pwkb only involves the negative frequency 
solutions of the Dirac equation. Thus this formula corresponds precisely to the naive 
picture of the Dirac sea as being composed of all negative-energy solutions of the Dirac 
equation. Theorem 15 . 1 1 and Theorem 15 . 21 show that the fermionic projector agrees with 
this naive picture, up to error terms which we now discuss. We first point out that, 
according to (]5.5p and (]5.6p . the fermionic projector has the naive scaling 



P(4>) ~ / <t>(f) R(t) dr . 
Jo 
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In order to compare with the error estimate (|5.15p . we need to assume that is 
supported away from the big bang and big crunch singularities, so that 

fir fir 

/ 0(r) R(r) dr ~ i? max / 0(r) dr . (5.16) 

This assumption is reasonable because we cannot expect the WKB approximation 
to hold near the singularities (in particular because "quantum oscillations" become 
relevant; see [14] ). Under this assumption, the estimate (|5.15p can be translated to 
a relative error of the order 0((mi? max )~5 ). We conclude that the error terms are 
under control provided that the size of the universe is much larger than the Compton 
scale 1/m. One should keep in mind that our theorems hold for a fixed function g 
in (|5.2p . This implies that the metric must be nearly constant on the Compton scale. 
Note that our estimates do not involve time integrals over the error, as one would get 
in a Gronwall estimate. This means that the local errors in different regions of space- 
time do not add up; we merely need to keep the error small at every space-time point. 
We also point out that, even when evaluating away from the singularities (see (|5.16p ). 
the behavior of the metric near the singularities still enters our construction via the 
integral (|5.10p . It is a main point of our analysis to estimate this integral without 
making any assumptions on the asymptotic form of g near the big bang or big crunch 
singularities. 

We finally discuss how our estimates depend on the momentum A. In view of (|5.14p 
and the error term in Theorem l5-H we may choose the quotient |A|/(mi? max ) arbitrarily 
large. This makes it possible to even describe ultrarelativistic Dirac particles. However, 
the constant c in ()5. 15[) and the error term in Theorem 15.11 depend on this quotient. 
This means that we cannot take the limit |A| — > oo for fixed mR max . It is not clear 
whether in this limit, the WKB approximation of P really breaks down or whether our 
estimates are simply not good enough to give a proper description of the corresponding 
asymptotic behavior. 

5.1. Computation of Swkb and P WKB . We now derive asymptotic formulas for 
Swkb and Pwkb including error estimates. 

Proposition 5.3. For any tq £ (0, it) there is a constant c which depends only on tq 
and the function g such that the matrix Swkb os defined by ()5.11j) has the explicit 
approximation 

\2 



SwKB = {[m iT ) V(T ° rl G -°i) r(To) + E 

with an error term E bounded by 



(5.17) 



c 



\\E\\ < — (5.18) 
m 

(here J.j is some norm on 2 x 2-matrices). Moreover, the eigenvalues [i^kb °f the 
matrix S W kb are given by 




f 



where 



T ):=-2 f \A 2 + m 2 R 2 . (5.20) 

J rn 
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Proof. A straightforward computation gives 

> N "T,t _ mR{T) 2 



WKBJ 



_l 1 ^WKB 

R(r) 



1 



R(t) 

+ -^rr4 COS . 



A 



f(r) /(r ) 
A mR(ro) 
-A 



mi?(r ) 
-A 



R i T ) ■ J. X 

+ -^H- sin <j> A 



mR(T ) J 

-i 

1 



Sir) """ r /(r ) Vm^o) "A y ' /(r) 

Carrying out the integral in (|5.11|1 . we can compute the eigenvalues of the resulting 
matrix to obtain (|5.19p . In order to derive asymptotic formulas, one must keep in mind 
that the factors sm<p and cos0 oscillate, resulting in small contributions to S*wkb- Let 
us quantify this effect for the integral involving cos<^> (for the integral involving sin</> 
the argument is exactly the same). We first transform the integral by 



cos 6 dr 



R(t) d 



— — sin <b dr 



d 



n 27(7) 2 dr 



— sin 6 dr . 



/o Sir) r J Q /(r)0'(r) dr 

Integrating by parts and using that R vanishes at both end points, we obtain 



io Sir) 
This yields the estimate 



cos 6 dr 



77 R (A 2 - m 2 R 2 ) 
(A 2 + m 2 i? 2 ) 2 



n R{r) 
fir) 



cos 6 dr 



< 



R 



A 2 + m 2 R 2 



dr 



\Xm\ 



sin 6 dr 



d ( mR 

— — arctan — — 
dr \ A 



dr 



On an interval where R is monotone, we can carry out the last integral, giving at 
most 7r/2. Since R is piecewise monotone, we can subdivide the interval (0, it) into iV 
subintervals on which R is monotone and carry out the integral on each such subin- 
terval. We conclude that 



R(r) 



cos 6 dr 



< 



1 Ntt 
|Am| ~ 



/o fir 

Next, a direct calculation shows that the matrix 

1 / A m 
7R) \mR(r ) 

has eigenvalues ±1 and is thus uniformly bounded. This completes the proof. 



R(ro)\ 



□ 



Proof of Theorem \5.1[ Writing the spectral calculus with residues, we have 



-X(-oo,o) 



WKB I 



j> (Swkb - A) 1 d\ 



where T is a contour which encloses the negative eigenvalue of 5 WK b- Estimating the 
integral in (|5.17p by 



mR 2 



o VA 2 + m 2 R 2 



dr > 



mR 2 



o ^A 2 + m 2 R 2 1 , 



dr = c- mR ™* 



^/X 2 + m 2 R 2 a , 



with 



c:= g z dr > , 
/o 
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we find that T can be chosen as a circle with center ju WKB and radius r given by 

mRl 



"-max 



?(o) 



Denoting the first summand in (|5.17p by S^kb and computing the contour integral 
gives 

-x(-oo,o) ft) = y^y 1 fg -1) F ( r °) • 

In order to estimate the error term E, we write the corresponding contour integrals as 
jf [(Swkb - A)" 1 - (S$ KB - A)- 1 ] dA = jf ^ + ^ - A)" 1 didA 



r J o 



^WKB 

Taking the absolute value and estimating the integrand, we obtain the error bound 



"max 



2ir Jy r 2 r m 2 R 2 n 

where in the last step we applied (|5.18p . Using ()5,12p . ()5. 13[) and (|5.9p gives the 
result. □ 

5.2. Estimates of U — £/wkb and S — S*wkb- The goal of this section is to derive 
the following estimate. 

Proposition 5.4. For any To 6 (0, it) there is a constant c which depends only on tq 
and the function g such that for all m and R max with mR max > 1 , 

_1 i 

\\S - S WKB \\ < cm 5it> max (5.21) 

(where \\.\\ again denotes a matrix norm). 

In preparation, we begin with three technical lemmas. Note that, as the matri- 
ces U T,T ° and C/wkb are both unitary, instead of JJ T,T ° — Uy^KB we can J us t wen estimate 
the matrix W(r) — 1, where W is the unitary matrix 

W(r):=(U^ B yU^. (5-22) 

A short calculation using (|5.3j) . (|5.7p and (|5.9p shows that 

d T W(r) = (U^ B y V(ry (d T V(r)) . 

Again using the definition of W(r), we obtain the differential equation 

d T W(T) = X{r)W(T) with X :=(U^ B )*V*(d T V)U^ B . (5.23) 

A straightforward computation gives 

_ Ami? l_ ( -«Asin(0) /o«>s(0) - imR sin(</>)\ 

2p f \-f Q cos^)-imR sm((j)) iAsin(0) )> [0 > 

where <f> is again the function (I5.20p and 

f := /(r ) , R := R(r ) . 
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Lemma 5.5. Assume that the function R(r) is monotone on the interval [ti,T2] C 
(0, 7r). T/ien 



||W(r) 



< 



1 



arctan 



mR(r) 
A 



Proof. Using Kato's inequality together with the fact that W is unitary, we know 
from (pT23l) that 



IW-II 



T2 



< 



Tl 



r-2 



\\d t W(r)\\dT< / \\X(r)\\dr. 

I Tl J Ti 

The matrix appearing on the right hand side of (15.24H is anti-Hermitian with eigen- 
values ±ifo- Hence 



\X\\ < 



XmR 




2/2 





1 d (mR\ 
— — arctan ( — ; — I 



2dr 



V A 



(5.25) 



where the last step is immediately verified by computing the derivative of the arctan 
and using (15.8p . Integrating on both sides and using that R is monotone gives the 
result. □ 



Lemma 5.6. For any tq G (0, tt) there is a constant c depending only on tq and the 
function g such that 

c\X\ 



1 < 



\\W(r) 

f'TT 

/ \\W(t) -1\\R(t) dr < 
Jo 



mR(r) 
cir\\\ 



m 



(5.26) 
(5.27) 



Proof The inequality ()5.27|) follows immediately by integrating ()5.26p . For the proof 
of (|5.26p . it suffices to consider the case r > To, because the case r < r is analogous. 
We choose intermediate points r\ , . . . , tn with 

T < Tl < • • • < TAT = TT , 

such that R restricted to the subintervals [Tfc_i,Tfe] is monotone for all k = 1,...N. 
Then r lies in one of the subintervals, r G [r n _i,T n ]. Applying Lemma 15.51 on the 
interval [tq,ti] and using that W(tq) = 1, we obtain 



2||W(r 1 ) - 1|| < 



fmR{r) 

arctan I — 

V A 



TT 



< I arctan , 

2 V A ) 



/ mR(T ) \ 



+ 



tt /mR(ri)\ 

arctan | . 

2 \ X J) 



Applying the elementary inequality 



gives 



TT 1 

arctan(x) < — 

2 x 



2||W(ri) - 1|| < 



for all x > 
|A| , |A| 



+ 



mR(To) mR(T\) 
Proceeding similarly on the other intervals, we conclude that 

|A| , , |A| ; |A| 



\\W(r) 



1 < 



17lR(T ) 



+ 



+ 



mi2(r, 



+ 



n-lj 



mR(r) 
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Using the scaling (|5.2p . we obtain 

||^(r)-l||<-^L (J_ + ...+ _J_ + 1 y 
mR{T) \g{T ) g{T N ) J 

giving the result. □ 

4 

Lemma 5.7. Suppose that A > (mi? max )s . T/ien i/iere is a constant c which depends 
only on tq and the function g such that for all m and i? max with mR max > 1 , 

\\W(r) -1 1| <c(miW)-B (5.28) 

r 71 " i i 

/ ||W - 1|| R(t) dr < cm-sR^ ax . (5.29) 
Jo 

Proof. We write (j5.24jl as 

X (r) = /i(r) A M (r) , (5.30) 

where 

Ami? 

h = and 
4/3 

1 ^ — zAcos(0) — /o sin(</>) — imRo cos( 



/o \/o sin(</>) — imi?o cos(0) iAcos(<^>) 
Integrating (|5.23|) . we can employ ()5.30p and integrate by parts to obtain 

T-2 



T 2 /• . . d 

ri 



(W(t) - 1) = / XW dr = hMW - J M jp(hW) dr 



n 

T2 



hMW 



n 



2 (jWW + MhXW ) dr . 

Tl 



where in the last line we used (|5.23p . The matrix M is anti-Hermitian and has the 
eigenvalues H. Moreover, the matrix X can be estimated by the first inequality 
in (|5.25p . which we now write as \\X\\ < \2fh\. Using furthermore that W is unitary, 
we obtain 

\\W(t) - 1||| T2 < \h(T 2 )\ + \h(n)\ + 2 + \2fh 2 \^jdT. (5.31) 

Now suppose that |A| > (m.i? max )^ with (3 < 1 (choosing (3 = 4/5 later will give the 
result). Using the estimate / > A, we obtain 

\h\ < R max \g\ < (m/W) 1 ^ \g\ 
\2fh 2 \<^RL x \g\ 2 <l(mR max r^\g\ 2 

\h\ < \ (mR^ x ) 1 -^ \g\ + | (mR m!lx f-^ \g\ 2 . 

Using these inequalities in (j5.31|) for t\ = To or t 2 = to, we conclude that there is a 
constant c as in the statement of the lemma such that 

\\W(r)-l\\ <c(m/W) 3 - 4/3 . 
Choosing (3 = 4/5 gives (|5.28p . Integrating yields (|5.29p . □ 
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Proof, of Proposition [531 We first derive a bound on the norm of S — .Swkb in terms of 
\\U — CAvkbII- Introducing the notation Y(r) = U T ' T ° (U^}^ b (t))* and applying (|5.10p 
and (|5.1ip . we find 



S-S- 



WKB 



1 
-1 



Y{r) 



1 




1 Y(r)* 



[7 T ' T ° R(t) dr. (5.32) 



Since U and Y are unitary matrices, (|5.32j) implies that 



\S-S- 



WKB 



< 



1 
-1 



r ( »(J \)Y(r) RAtXIt 



< 



£\\(1-Y(r))(l ^~Y(r)(l (1 - Y(r)*) || R(r) dr 

(||l-Y(r)|| + 111 -Y(t)* II) R{r)dr 



o 



2 / £T™ - E^(r) J2(t) dr = 2 / || IU(r) - 1 i2(r) . 
'o jo 



Now, Lemma 15.61 yields (|5.2ip . while Lemma 15.71 gives the remaining case. This com- 
pletes the proof. □ 



5.3. An Estimate of P — Pwkb- 

Proof of Theorem \5.$\ Introducing the abbreviations 

■A/" = —X(-<x>,o)(S) an d A/Wkb = —X(-<x>,o)(Swkb) , 
we obtain from (15. 5h and (I5.12p 

P — P\VKB = -A/" k m — A/\VKB &WKB 

= (Af — A/"wkb) k m + A/\vkb (km — k WKB ) . 

Applying a test function cp £ Cq°((0, 7r), C 2 ) and taking the norm, we can use that A/"wkb 
has norm at most one to obtain 



{P-P^ 



WKB 



< 



in — "-WKB 



+ ||A/"- A/wkbII 



(5.33) 



In order to estimate the first summand in (]5.33p . we first note that, according to (15.6 
and (15331) . 



(Ar m - fc W KB)(0) 



2tt 



(tt t ' t o _ rr 1 "' 7 " "i* / 

1 ' u wkbJ \ Q _2 



(r) fl(r) dr . 



(5.34) 



Using (15.221) and Lemma EH we get the estimate 

m J 



\{k m - fcwKB)(0)|| < 



T)\\dT. 



In the case |A| < (mi? max )s, we get the inequality 



\(k m - k WKB )((j))\\ <c(mR 



max J 5 ^''ir.ax 



(r)||dr. 



(5.35) 



In the remaining case |A| < (mR max )5, we apply Lemma [5.71 to (|5,34p . again giv- 
ing (^351) . 
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It remains to estimate the second summand in f|5.33j) . Noting that 

HM0)ll < ^ f U(r)\\ R(r) dr<±- J* ||0(r)|| dr , 
the proof of the theorem is completed by applying Lemma 15.81 below. □ 



Lemma 5.8. Under the assumptions of Theorem \5.S\ 

\W -MwkbW < c(mfi max )"5 . 
Proof. Writing the spectral calculus with residues, we have 

2m J T 

where T is a curve in the left half plane enclosing all negative eigenvalues. Similarly, 

A/"wKB = 77— / (SwKB ~ 0^ d( ■ 

We choose T WK b as a circle centered at the negative eigenvalue // WK b with radius 
r = |a* WK bI/2- Using (|5.19p together with (I5.14p . we can estimate this eigenvalue by 

r mR 2 

l/wl -/ twt^w 



mR? 

o y / fc 2 m 2 i?^ iax + m 2 R 2 



> I , _ _ _ dr = cR max , 



c ■- I , g2 dr > . 



where 



Vk 2 + g 2 

According to Proposition 15.41 we can treat the operator AS := 5 WK b — S as a 
perturbation. More precisely, the min-max-principle (see for example [23]) yields that 
the negative eigenvalue of the operator S + tAS, t £ [0, 1], lies inside T, and that the 
distance of the eigenvalues of all these operators from T is at least equal to 

d := - > . (5.36) 

2 ~ 4 v ' 



It follows that 

AA-A/Wkb = 

2m 



(s - cr 1 - (s WKB - cr 1 
^l(j t {s+tAS - crldtdC 

— if (S + tAS - Cr 1 AS (S + t AS - C)" 1 dt d( , (5.37) 
2m Jr Jo 

where we set AS = SVkb — S. Taking the norm and estimating gives 

||AA-A/wkb|| < max / 11(5 + tAS - O'H 2 ||A5|| d|C| < 4 ||A5|| . 
27Ttg[o,i]j r « 

Applying (|5.36p and Proposition 15.41 gives the result. □ 
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6. Discussion of Klein's Paradox, an Example where 5 = 

Qualitatively speaking, the results of Section [5] show that our definition of the 
fermionic projector reduces to the naive notion of the Dirac sea as "all solutions 
of negative frequency," provided that the metric is nearly constant on the Comp- 
ton scale. This raises the question what happens if the metric varies substantially on 
the Compton scale. In order to tackle this question, we now analyze the situation 
for a closed Friedmann-Robertson- Walker space-time with a piecewise constant scale 
function. This analysis is also instructive because it will give a connection to the 
well-known Klein paradox (see for example Section 3.3] or [25j Section 4.5]). 

We again consider the line element (|5.1[) . Again separating the spatial dependence, 
the operator S is given by (|5.10|) . where the unitary matrix U T ' T ° is defined as the 
solution of the initial value problem (|5.3p and (|5.4p . In order to get a better geometric 
understanding of the dynamics, it is useful to decompose the matrix in the integrand 
of ()5.10p in terms of Pauli matrices by setting 

1 m I ttQ,t f 1 \ rrr n 



v a (r):=-Tr\a a U^\: , ) • 

After cyclically commuting the factors in the trace, we obtain 

v a (r) = (w a (r),e 3 ) , (6.1) 
where uT a (r) is (for any given a = 1,2, 3) the vector 

tf a (T) = ±Tr(BU%°<T a U% r ). (6.2) 



Taking the r-derivative and using ()5.3[) gives 



d T w a (T) = ^Tr ( d 



^-Tr([a,da]U^°a a U, 



0.T 

m 



2 

where the vector d has the components 

d = 2 (A, 0, —mR) . (6.3) 

Using the commutation relations of the Pauli matrices, we obtain 

d T w a = d A w a . (6-4) 

Moreover, evaluating (|6.2p for r = gives the initial condition 

w a (0) = e a . (6.5) 

The differential equation (|6.4|) describes a rotation of the vector w around the axis d, 
which also depends on r. This equation can be regarded as the Block representation 
of the Dirac equation (|5.3p (see the discussion of the Dirac equation in p~4j Section 2]). 
However, the initial conditions (j6.5j) and the connection to the vector v by (|6.ip . are 
specific to the construction of the fermionic projector. 

Next, we choose the scale function R(t) to be piecewise constant. Thus we introduce 
intermediate points tq = < t\ < ■ ■ ■ < tn = vr and set 

N 



71=1 
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with parameters R\, . . . , Rn > 0. Then on the subinterval [r n _i,T n ), the dynamics 
of (|6.4p reduces to the rotation of the Bloch vector around the fixed rotation axis 

d n : = 2 (A, 0, —mRn) . 

The angular velocity of this rotation is given by 2\J A 2 + m 2 i? 2 . This is the frequency 
of the so-called Zitterbewegung of the Dirac particle; it is twice the frequency of the 
oscillations of the Dirac wave functions. We denote the number of full rotations of the 
Bloch vector on the interval [r n _i,r n ) by p n . Then 

_ _ TTPn 

Tn Tn ~ l ~ ^/X 2 + m 2 Rl ' 

If the scale function is constant, we may decompose the spinors into eigenfunctions 
of the matrix da. This corresponds precisely to the splitting of the solutions into 
solutions of positive and negative frequency. However, this splitting depends on the 
value of the scale function. In particular, if R changes discontinuously, the canonical 
splitting into positive and negative frequency solutions gets lost. This is very similar 
the situation discussed in Klein's paradox (with the only difference that we here con- 
sider a gravitational instead of an electromagnetic field; also usually a potential with a 
spatial dependence instead of a time dependence is considered). Let us analyze what 
happens in the construction of the fermionic projector. According to (|5.10p . we must 
integrate v over time, 

fir 

S= v{j)R{r)dr. (6.6) 
Jo 

As a consequence, only the "time average" of v enters the construction, but a canonical 
splitting of the solution space into solutions of positive and negative frequency is no 
longer needed. In the situation of Klein's paradox, this time average means that the 
fermionic projector will be an "interpolation" of the concepts of negative frequency 
before and after the step potential (for a related discussion of a scattering process 
see Section 5]). This interpolation is performed in such a way that the construction 
of the fermionic projector is manifestly covariant and independent of observers. 

According to (|6.4j) and (j6.3|) . the Bloch vector w rotates around a time-dependent 
rotation axis d. This can lead to bizarre effects when the rotation axis is tilted several 
times in a fine-tuned way. In order to illustrate such effects in a simple example, we 
now construct a space-time where 5 = 0. In this case, the fermionic projector defined 
by (15. 5p vanishes identically. By slightly changing the geometry, one can perturb the 
eigenvalues of S. The operator — X(-oo,o)(>>) i n fl5.5|) is certainly not stable under 
such perturbations. This shows that in a space-time with S = 0, the definition of 
the fermionic projector suffers from an instability and thus depends sensitively on the 
detailed geometry of space-time. From the physical point of view, this shortcoming 
does not seem to be of any significance, because the class of space-times with 5 = 
seems quite pathological. 

We now introduce our example in detail. With the scale function, we can adjust the 
angle of the rotation axis d to the z-axis. We choose i?i and i?2 such that this angle 
equals 10° resp. 70°, i.e. 

R 1 = — cot(10°) , R 2 = — cot(70°) . (6.7) 
m m 

Moreover, we always choose the parameters p n such that we rotate a half-integer times 

around d n , so that the rotation amounts to a reflection at the axis d n . Composing 
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Figure 1. The functions v a (r). 

the reflection at d\ with the reflection at ofe gives rise to a rotation around the axis e2 
about an angle of 2 • 60° = 120°. Repeating this construction three times gives in total 
a rotation by 360°. More specifically, for the construction so far, we choose N = 6 and 

(Rn)n=l,...,6 = (Ri,R2, R u Ba,Rx,Ri) , (Pn)„=l,...,6 = (5-5, 0.5, 5.5, 0.5, 5.5, 0.5) 

with R\ and R2 as in fj6.7j) (the choices p\ = 5.5 and P2 = 0.5 are arbitrary; other half- 
integer values would work just as well). Moreover, for convenience we chose A = 3/2 
and m = 1. We solve the system (|5.3|) and (|5.4j) with to = 0. Decomposing the 
resulting signature operator S, ()5.10|) . in terms of Pauli matrices (|6.6p . it follows by 
symmetry in the ei/e3-plane that S\ = S3 = 0. However, S2 7^ 0. 

In order to also arrange that S% = 0, we take the above space-time twice, with the 
opposite time orientation. Thus we now choose N = 12 and 

{R n )n=i,...,i2 = (Ri j R2 ) Ri , R2 , R\ , Ri , R2,Ri,R2,Ri,R2,Ri) 
(p n )n=i,...,i2 = (5.5,0.5,5.5,0.5,5.5,0.5, 0.5,5.5,0.5,5.5,0.5,5.5) . 

Then the symmetries imply that 5 = 0. To illustrate the construction, in Figure [1] the 
functions v a {r) are plotted. 

Acknowledgments: We would like to thank Olaf Miiller and Jan-Hendrik Treude for 
helpful discussions. 
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